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Coefficient inaccuracy in transversal filters is known to degrade 
the frequency response, particularly in stopband regions. The mag- 
nitude of the stopband degradation increases with the number of 
stages n, the length of the impulse response. A widely used for- 
mula for the error in frequency response is proportional to vn. By 
extending recent results on random trigonometric polynomials, we 
show that for random additive coefficient errors with variance a 2 , the 
error hH(a) in frequency response for large n is such that 



max | &H(o}) | « oJn log n 

where log denotes the natural logarithm. This result leads to an 
absolute bound on attainable stopband rejection for any band-select 
transversal filter with given coefficient inaccuracy. In particular, the 
result places a definite limitation on the quality of band- select filter- 
ing that can be achieved using a ccd split-electrode filter. It also 
implies bounds for the peak sidelobes of random radar arrays. 



I. INTRODUCTION 

In recent years, the transversal filter has emerged as an essential 
signal-processing structure for a large variety of applications in com- 
munication systems. A few of these applications are matched filtering 
in radar or spread-spectrum systems, equalization in data receivers, 
echo cancellation for satellite communications, and band-select digi- 
tal filters. The term "transversal filter" originally referred to the 
continuous-time tapped delay line structure where an output is formed 
from a weighted sum of the tap voltages. The same basic function has 
also been achieved using lumped networks to approximate the delay 
sections. More recently, transversal filters have been realized with 
digital circuitry using shift registers and digital multipliers, operating 
on a sampled and quantized input signal. The most recent development 
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is the emergence of two new technologies, charge-coupled devices 
(ccds) and surface acoustic wave (saw) devices which allow the 
realization of discrete-time transversal filters without the need for 
analog-to-digital conversion. 

The new technological advances now offer the possibility of realizing 
transversal filters with hundreds and perhaps even thousands of tap- 
weight stages on a single integrated-circuit chip. These developments 
suggest that extremely sophisticated signal-processing functions can 
readily be obtained. Specifically, with a sufficient number of taps, a 
transversal filter can be designed to approximate virtually any specified 
frequency response as closely as desired. However, the inevitable 
inaccuracies in implementing the desired weighting coefficients result 
in a departure of the actual frequency response from the predesigned 
frequency response which increases with the number of tap-weight 
stages. In digital filtering, coefficient values can be made as accurate 
as needed, but at the price of increasing hardware costs. With the ccd 
or saw technologies, there are fundamental limits on attainable accu- 
racy. Also, in adaptive filtering, the weight-adjustment algorithm 
results in a steady-state coefficient inaccuracy. It is therefore necessary 
to have a quantitative knowledge of the degradation in performance of 
the transversal filter as a function of the coefficient inaccuracy and 
the number of stages. 

For most applications, the appropriate performance measure for the 
realized transversal filter is the maximum deviation in frequency 
response magnitude from the desired values over the particular fre- 
quency band of interest. In this paper, we focus on this performance 
measure by examining the error-frequency response due to coefficient 
inaccuracy and show that under very ge neral co nditions the maximum 
magnitude is given asymptotically byaVn log n, where n is the number 
of stages, a is the rms coefficient inaccuracy, and log denotes the 
natural logarithm. Several other closely related results and implica- 
tions are also presented. 

Since the attainable quality of a designed filter increases with n, the 
number of stages, and for a given coefficient inaccuracy the degrada- 
tion increases with n, the question arises: Is it possible to realize a 
filter with arbitrarily high quality in spite of a given coefficient inac- 
curacy if n is made sufficiently large? We make this question more 
precise later and show that the answer is negative for low-pass filtering 
with a transversal filter structure when "quality" is measured by the 
amount of stopband rejection. In other words, a limit on filter accuracy 
implies a limit on attainable filter quality regardless of the number of 
stages used. The results of this paper provide a tool for determining 
the ultimate limitation on transversal filter performance associated 
with a particular technology or a particular adaptive algorithm for 
weight adjustment. 
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In ccd transversal filters, the split-electrode method requires that 
the tap weights be scaled so that the maximum magnitude of the 
coefficient values is unity. The pattern generator used in making the 
photomasks for ccd fabrication introduces a quantization error whose 
peak size is a fixed fraction of the maximum coefficient magnitude. 
Now, for most applications, increasing the number of stages to be 
realized corresponds to including additional coefficient values repre- 
senting the tail of the desired impulse response. Consequently, increas- 
ing n does not alter the scaling of the coefficient values for ccd 
implementation. As a result, a coefficient error can indeed be modeled 
as an additive random variable whose variance does not depend on the 
desired coefficient value. 

A problem that is very similar to that considered above occurs in 
the theory of random arrays. 1 These are arrays consisting of fewer 
elements than conventional phased arrays, with the locations of the 
elements in the array picked randomly. Such arrays are less costly 
than conventional phased arrays, but this advantage is gained at the 
cost of increasing the peak sidelobes. Our main result shows how big 
those sidelobes can be expected to become. 

II. PROBLEM FORMULATION 

Regardless of the particular application, the transversal filter may 
be described by its frequency response, //(to), which has the general 
form 

//(to) = 2 a***, (1) 

where to is t he n ormalized frequency variable, L is the number of 
stages, j = v— 1, and the coefficients a* are real-valued numbers 
specified by the designer. Since //(to) is periodic, only the frequency 
interval < to < 2tt- need be considered. 

We note, in passing, that (1) also describes the discrete Fourier 
transform, so that the results of this paper are also applicable to 
studying the effect of approximate representations of given data values 
on the Fourier transform of the data. 

A special case of transversal filters, of particular interest in band- 
select filter design, arises when the coefficients are chosen to have the 
symmetry property: 

a A -«L-*-i for 0S*<(I-1). (2) 

When L is odd, this condition results in a linear phase transfer function 
having the form 

n 

//(to) = e +Jun £ b„ cos ku, (3) 
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with n = (L - l)/2, and 

bk = 2 a„-* for k j^ 0, bo = a„. 

Implementation of the coefficients a* for the general form (1) or bk 
for the linear phase form (2) inevitably results in the introduction of 
errors or inaccuracies. We denote the actual (inaccurate) value realized 
as a'k, or as b% for the linear phase case. Then the Ath coefficient error 
is the difference e* = a* — a*, or e* = b'k — bk in the linear phase case. 
The realized transfer function then differs from the desired transfer 
function by the error transfer function defined as 

AM = 2 ***- (4) 

in the general case or, in the linear phase case: 

g n (u) = e +Jun i €k cos ku>. (5) 

Clearly, the error transfer function, if known, provides a full description 
of the degradation in performance of the realized filter from the desired 
performance in the absence of inaccuracies. 

Since the errors, €*, are generally not known prior to fabrication of 
the filter, they are modeled most effectively as random variables whose 
distribution depends on the particular mechanism involved in fabri- 
cating the tap weights. In digital filtering, the errors are due to 
coefficient quantization and are usually modeled as uniformly distrib- 
uted random variables. The error terms for different coefficients, being 
independently produced, can reasonably be assumed to be independent 
random variables. 

Additive error components were used by Knowles and Olcayto 2 for 
modeling coefficient quantization in recursive filters. Chan and Rabi- 
ner 3 applied this approach for transversal filters and evaluated the rms 
values of /l(w) and g„(u) at a particular frequency. They assumed 
mutually independent and uniformly distributed errors €* resulting in 
rms values for the error transfer function proportional to vX, or Vn in 
the linear phase case. By taking the maximum over all frequencies of 
the rms deviation, a frequency-independent upper bound on the error 
transfer function is obtained which is valid at any particular frequency 
with high probability. 

More recently, Heute 4,5 noted that the bounds of Chan and Rabiner 
underestimate the degradation due to the maximum of \g„(u) | over 
the frequency band. It is this latter measure of degradation that is 
meaningful in most applications. Chan and Rabiner's bound is not a 
high probability upper bound for the maximum ripple magnitude taken 
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on by gn(u)- Heute proposed a heuristic upper bound for the maximum 
of | g n M | which has the form Q[a + bn + (en + d) I/2 ], where a, b, c, 
and d are constants and Q is the peak amplitude of the uniformly 
distributed error terms e*. His bound gave an improved fit to simulated 
data for values of n up to 128. We shall see later that Heute's bound, 
which for large n grows linearly with number of stages n, grossly 
overestimates the degradation as n becomes much larger than 100. 

Andrisano and Calandrino 6 assumed that the error transfer function 
is a Gaussian process and found an (implicit) bound on stopband 
rejection as the solution of a transcendental equation. 

In this paper, we take as the measure of degradation due to coeffi- 
cient inaccuracy, 

DL = max\f L (u)\ (6) 

ucQ 

for the general transversal filter and 

M n = max \g„(u)\ (7) 

ucQ 

for the linear phase transversal filter, where ft is a particular frequency 
band of interest. We assume the errors e k are mutually independent 
random variables with a common distribution satisfying certain regu- 
larity conditions that include the uniform and normal distributions as 
special cases. 

We establish here for the first time that the maximum frequency 
respo nse erro rs D„ and M„ are asymptotically (for large n) given 
by oVn log n where o is the rms coefficient error. Although the result 
is asymptotic, Lawrence and Salazar 7 found that it was moderately 
accurate in one study of a low-pass filter with only 33 taps. Application 
of the result to low-pass filter performance is examined briefly in this 
paper and more extensively in Ref. 8 and 9. Until the report of our 
result," the correct behavior of the error frequency response magnitude 
had apparently not been recognized in the digital filtering literature. 

The existing mathematical results most closely related to our work 
are due to Halasz'" who considered random trigonometric sums with 
coefficients that take on the values ±1 with equal probability. While 
too restrictive to apply to transversal filters, his methodology was 
useful in deriving our more general upper bound on the maximum 
error frequency response. 

Our main result is also applicable to the analysis of random arrays, 
and in particular to that of statistical arrays. 1 These are arrays con- 
sisting of k isotropic radiators placed among n positions (n > k) that 
are spaced A/2 apart (A = wavelength), with the k positions to be 
occupied by the k elements determined at random. The array factor of 
such an array is defined as 
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flu)-! *■**". (8) 

r-0 

where g r = 1 if the rth position is occupied by a radiator, and g r = 
otherwise. This can be rewritten as 

f(u)=- J rf*"+J ere^", (9) 

1 r-0 r-0 

where € r = 1 — £/n for the A values of r for which g r = 1, and e r = 
— k/n otherwise. The first sum above represents (except for the k/n 
multiplier) the array factor of a conventional phased array. The 
random choice of the positions for the radiators corresponds to letting 
the e r be independent random variables, assuming the value 1 — k/n 
with probability k/n, and the value — k/n with probability 1 — k/n. If 
we assume that k ~ an as n — > oo, then our theorem shows that, with 
the probability approaching 1 as n — > o o, the se cond sum in (9) will 
never be significantly larger than Vl— a Vn log n and that, conversely, 
it will get that large on any subinterval. This result, which had been 
derived only heuristically before, 1 explains why random arrays are 
usually not very satisfactory. 

III. STATEMENT OF MATHEMATICAL RESULTS 

As we saw in Section II, the errors in the realized transfer functions 
are given by £*-o €* e Jkw , or £*_o e* cos ku, or £*-o e* sin kw. Hence the 
distribution of the random variables €* will depend on the model for 
the sources of inaccuracies. For digital implementation, the usual 
model assumes e* to be independently distributed uniformly between 
—A, +A, where A is the maximum error due to truncation of the 
coefficients of the filter. In other situations, a Gaussian distribution 
may be more appropriate. But, as we shall see, the asymptotic behavior 
of the maximum magnitude of the error is not dependent on the exact 
nature of the distribution. It depends only on a few functional of the 
distribution. 

The results presented here rely on an important assumption about 
the distribution of the €*. We assume throughout that the e* have 
mean zero and finite sixth moment, so that the characteristic function 
E(e ixt ") of e k is such that 



£(*/**) = exp 



- X arf + O (x 6 ) 



(10) 



for x in some nontrivial interval [-d, d]. (Note that a 2 > if the e* are 
not identically zero.) Condition (10) is satisfied for most probability 
density functions of practical interest. 
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Now we are ready to state the main result: 

Theorem: Let c*, k = 1, 2, • • • be a sequence of independent identically 
distributed random variables satisfying (10). 



Then there exist constants C\ and C>, not depending on n, such that 



max 



< -J2a> sin log n + C\\ log log n 

V log n 



holds with probability > 1 — C> (log n)~ A . Furthermore, if £2 is any 
subinterval of[0, 2tt\ of length > (log n)~ { and a is any real number, 
then 



max 

dtV. 



Re I e j " £ e k e kjH \ > J2a~> sin log n - C, Wp- log log n 



holds with probability > 1 — d (log n)~ 4 . 

Thus, with high probability, max | f(9) | is about \l2a-i sin log n. The 
proof is outlined in Section V. 

Remark 1. By choosing a appropriately, we can conclude that each of 
£ €k cos (kd), X e* sin (kd) becomes large on any long 6 interval with 
high probability. 

Remark 2: The estimates presented here are not the best possible 
ones. For example, the interval £2 in the lower bound proof can be of 
size n (log n)~ x for any A > 0. 

IV. APPLICATION TO LOW-PASS FILTERS 

The usual specifications for fir low-pass filters are shown in Fig. 1." 
A design problem is to find the smallest n such that 

71-1 

^] Ok cos kd 

lies between 1 — 5i and 1 + 8\ in the passband, i.e., for e [0, F p ] and 
between and 8> in the stopband, i.e., for e [F s , ir\. Estimates for n 
given Si, 5 2 , F p , F s are given in Ref. 12. However, the validity of the 
estimates in Ref. 12 for regions of practical interest is not proven. An 
empirical relationship is given in Ref. 11. 

As n increases, smaller 5 2 , 8i and F s — F p are possible. Hence, a 
question that is usually raised is: Given that the a/,'s cannot be realized 
exactly, what can be said about the minimum 82 possible if the 
distribution of €*, the error in ot k , is known. If a*'s could be realized 
exactly, arbitrarily small values of 8> can be obtained by making n 
large. However, £//s introduce errors that grow with n as seen from the 
theorem. So there is a trade-off between errors introduced by inaccur- 
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Fig. 1— Frequency as fraction of 2 w. 



acies in a*'s and the improvement in performance with increasing n. 
For an example, we consider the stopband rejection, 20 logio 8 2 , as a 
figure of merit with given values of 6\ and F p , F s . The empirical 
formula gives the following relation for n, the minimum n required to 
achieve a stopband rejection of 20 logio 8 2 . 



n = d log 8 2 + c 2 , 



(11) 



where d and c 2 are constants depending on 8\ and F a — F p . u For fixed 
point digital implementations, if the coefficients of the filter, the a*'s, 
are truncated to d bits, then the "error" in a* is generally modeled as 
a uniform random variable €* having values between — 2~ d and 
2~ d = A. For this model, a 2 of the theorem is 1/6 A 2 . Hence, the 
maximum error e„ due to these inaccuracies, 



e„ = max 



n-1 

£ €k COS kO 

k=0 



is such that 



y/nlogn > 3 V3 



as n 



(12) 



(13) 
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and 

A 



„ = -Jn log n 

n/3 




<c \l- log log n (14) 

logn 



with probability > 1 - 0((log n)~ 4 ). 

Using the limit (13) to indicate expected deterioration in perform- 
ance, we can arrive at a design rule. If coefficients are truncated to 
d bits, then the minimum achievable 8 2 before the random errors 
become comparable to 62 itself is given by: 

2-rf 

—j= V(ci log 5 2 + c 2 ) log (ci log 62 + c 2 ) = 8 2 . (15) 

V3 



Putting 8 2 = 2 _s , 

2~ d 



Va V(— c\s + c 2 ) log (— c\s + c 2 ) 



(16) 



where c\ = C\ log 2. 



From the above formula, we can estimate the required precision for 
the coefficients for a given value of 62 = 2~\ 

In design of ccd filters, a similar formula can be used. In situations 
where the tap-weight errors can be modeled by a Gaussian random 
variable with a standard deviation A, then a 2 for our theorem is A 2 /2. 
The minimum achievable S 2 satisfies 

/ 82 

V(ci log 8> + C2) log (ci log 82 + c 2 ) = — . (17) 

Solving for 8>, we can estimate the optimum value of n. 

As an illustration of the effect of coefficient inaccuracy on limiting 
the stopband rejection of a low- pass filter, Fig. 2 shows how the best 
achievable rejection depends on the number of stages, n, for various 
values of the transition width AF = Fs — Fp. These curves were 
calculated by solving the empirical formula of Ref. 11 for 82 and adding 
to it the maximum error oVn log n. This gives an expression for the 
best attainable stopband rejection in the presence of coefficient errors, 
as a function of n, 8\, and AF. For additional curves obtained in this 
way, see Ref. 8. Computation also shows that varying the allowed 
passband ripple 5i has a negligible effect on the maximum attainable 
stopband rejection. It is evident that coefficient inaccuracy places an 
ultimate limitation on the attainable quality of a low-pass filter imple- 
mented with the transversal structure. 
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Fig. 2 — Best obtainable stopband attenuation for a low-pass transversal filter in the 
presence of coefficient inaccuracies. Root-mean-square coefficient error = 0.001, 
passband ripple = 0.0122. Curves are shown for four values of the transition width. Note 
that, for each curve, an optimum number of taps exist. Reducing the passband ripple 
allowable has the effect of shifting these curves to the right while reducing the peak 
value of each curve. 
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APPENDIX 

Here we outline the main steps in the proof, which follows that of 
Halasz, 10 in which he assumed e* to be ±1. (An earlier proof of a 
slightly weaker result had been outlined by Whittle. 13 ) Results that 
are incidental to the main line of reasoning are collected at the end of 
this outline. Let 

f(a, d) = I €„ cos (kB + a) = Re( e J " % €„e Jke 

(i) We construct a nonnegative function u(x) < 1 which can be 
used to indicate in an approximate sense the set of values of x that 
exceed given values. Let Mi, Mi, D > be given numbers. Then 
u(x) is zero for — M2 < x < Mi, and u(x) = 1 for x > Mi + D or 
x < -M% - D. In the interval [-M 2 -D, -M 2 ] and [M u Mi + D], 
u (x) is 40 times differentiate and u lr) (x) = 0(D~ r ) as D f 00, for < r 
<40. 
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For deriving the upper bound, we pr oceed a s follows: 
(i i) Put M i = Mo = M = V2a2 Vn log n + gD log log n where 
D = Jn/log n and g = 20^^72. 

(Hi) Let i;,(0 = 1/2tt /-. (1 - w(x))c- >te rfx. Then | m(t) \ = 0(M) 
and | t r v x (t) | = 0(D r+l ), 1 < r < 40. 

(iu) Let G = Jl" rfa /§" d6u(f(a, 6)) and u(f) = 8(t) - Vx{t), where 
5(<) is the Dirac delta function. Then 



G= da dd \ exp 



jt X e* cos (A0 + a) 



v{t)dt. 



Using the properties of v u we can then show (for details, see the 
section at the end of this proof) that 

E(G) = 0(n- i (\ognr U) ). 

{v) Let T = max | f(a, 6) \ . Then using the inequalities 

\2-f(a,$)\<*Tn, |f /-(a,0)|<r 
do da 

we can show that G < l/(n log~n) ==> T < M + 2D for large enough n. 
(vi) The result from step (iv) implies 

M Ga ^k) =o (o^)- 

Hence, using step {v) 

max I f(a, 6) I < >/2a^ V/i log n + (g + 2) \ log log n 

a.e V log n 

with probability > 1 - °( (1 ^ j • 

The derivation of the lower bound is more difficult, but similar. We 
will only outline the proof. We examine the values of f(a, 6) at the 
points 6 m = 2tt- (2m - l)/2n, for 1 < m < n. 

(vii) Let Mi = M = V2a^ sfnlogn - #Z) log log /i, and M 2 = 2M. 
Let S be a subset of the integers from 1 to n with cardinality greater 
than n(log n)' x and put F = ^ ( s«(/(ff, 0m)). 

As in the derivation of the upper bound, we can find the asymptotic 
behavior of the first two moments of F using the properties of u. 

(viii) We can show E(F) > c 3 |S|rc"'(log n)' n/2 for some constant 
c, > 0, and E(F 2 ) - E 2 (F) = 0(|S|n-'(log n) 21/2 ) + 0( \S\ 2 >r 2 (log 
n) 7 ). 

(ix) Now 

Pr{max f(a, d„>) > M or min /"(a, m ) < -2M} > 1 - Pr{F = 0} 

meS mtS 

COEFFICIENT INACCURACY IN FILTERS 2311 



by the definition of n. But Pr{min f(a, d m ) < -2M} = 0((log n)~ B ) 
from the upper bound, so that 

Pr{max f(a, m ) > M) > 1 - Pr(F = 0) - 0((log n)~ s ). 

mtS 

Further, Pr{F = 0} < Pr{(F - E(F)) 2 > E 2 (F)}, so by Chebyshev's 
inequality 

E((F-E(F)) 2 ) 



Pr{F = 0} < 



E 2 (F) 



(x) Using the bounds from step (viii), Pr{F = 0} = 0((log n)~ 9/2 ). 
Therefore, using the definition of Mi, M 2 , and D, we have: 



PrJ max f(a, 9 m ) > V2a 2 -Jn log n - g \ log log n \ 

[ m t s V log n J 

> 1 - Pr{F = 0} - 0((log n)~ s ) 

>l-0((logn)- 9/2 ). 

Details of Step (iv) 

From the definition of v(t) [see step (iv)], we can show that 



t\ r \v(t)\dt = 0(D~ r ) l<r<18 



J\l\>d/2 



t\ r \v(t)\dt = 0(D~ 19 ). 



Since the Fourier transform of t r v(t) isj r w r) (x), 



(18) 
(19) 



e-' ul t r v(t)dt\ = —^= | I u ir) (x)e- x2/4p dx 



= oU- x/2 D- r \ 



*'d» 



= o 



\MD r€ ) 



uniformly in f$ > 0, < r < 18. 
Similarly, for /? > 



J expf- 



•fi? I cos 2 (kd + a) ) t r v(t)dt = O ( -^ exp(-M 2 /Q) ) , 



[md 2 



(20) 



(21) 
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where 



Further, 



and 



Q = 40 £ cos 2 (kO + a). 



V cos 2 (Atf + a) = £ + - £ cos 2(*0 + a) 



£ cos2(^ + a)< T -r--^ T . 
k -i |sin0| 



*-i 



So 



n V0,a 



X cos 2 (*0 + a) < J ra „ 

*-i 1 - + 



fo *!°£IU,*|, 
2 2 log n 2 n 



since 



sin0 >r for - — — < 0|<*r-- — — . 

log n 2n In 



Therefore 



1 



-M 2 



exp 



d0 = 0(e- M2/2pn ), 



o 40 X cos 2 (kB + a) 

*-i 



whence 



J>J>£ 



exp 



-P? £ cos 2 (A:0 + a) \t r v(t)dt 



a) ful 



(22) 



n(^le- M2/2pn \ 
U \MD r€ )' 



(23) 



Using (10) and (23) above, we can derive the result of step (iv) as 
follows: 

Since €* are independent, 



£(exp 



jt X e * c * 

A-l 



= J] Ee JUkCk 
k-i 

-Y.attf i c' k + 0(n* 6 )| 

1-2 A-l J 
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= exp 



for \t\<d 



from assumption (10), where c k denotes cos (kd + a). Further, 



El exp 



jt £ € k Ck 



= exp 



- a 2 f s d 



+ £ ajt J £ c{ exp 
/-a k=i 



- a 2 t 2 I cl 



+ lliajt J l,ei\'vxp 



- a 2 t 2 X d 



*-i 



+ 0(nr) + 0(n 3 U| 9 ) for |*|<d, (24) 



since 



e - = e - b + (b- a)e~ b + ¥2 (b - a) V + 0( | b - a \ 3 ) 
uniformly for b e @, b > 0, a e % Re (a) > 0. Now we consider 



E 



exp 



jt £ e k c k 



r(!)di ] = J #(exp 



jt 2 € A c* 



!;(*)&; 



the expression on the right-hand side of (24) can be substituted for the 
integrand in the interval 1 1 | < d. Outside this interval, we can use the 
simple bound | exp [jt £? e k c k ] | < 1, and arrive at: 



E exp 



jt X € k C k 



v(t)dt 



-jV-p 



jt £ €kCk 



+ 0[ \v(t)\dt 

'I'M 



v(t)dt 



= I ffezp /t£«c* \)v(t)dt 
+ o(\ \v(t)\dt) 

\J\t\zd ' 

+ o(n I |*| 5 |i;(*)|tft) 
+ o(n 2 J t xo \v(t)\dt 
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+ 0(/i t 6 \v(t)\dt 



+ o(n 3 [ \t\ 9 \v(t)\dt\. (25) 

From (18) and (19), we see that the right-hand side is 

■Li 

-f 



jt y dkCk 



Ey exp 



jt Y e k Ck 



v{t)dt + 0(nD~ 6 ) + 0{n 2 D- }9 ) + o(n 3 D" 9 ) 

(26) 



vmt+ o{^4 



\ n 



To find the asymptotic behavior of E(G), we use (21). After integrating 
with respect to a,0, we have, for each of the terms in (24), with 
expressions in square brackets corresponding to those in (24), 



da I dO\ exp| - »/l a 

o Jo J-" 






M 



da\ dd £ aj £ c J k t J exp[ ]v(t)dt = Ol n—-. , 

Jo Jo y " 3 ' J-OD V 



-M i /2a i n 



da I dd \ (2 C/^Id. } exp| jrO'W,' 
o Jo 



J-3 



U \ n MTf e )' 



Therefore, collecting the previous results, we have 
E(G)=\ da\ dd fi(expl/*£c*CA 



v{t)dt 



. (^ e -«'^ + <^p) . 0(n -. (log „)-•). 
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